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A relativistic version of cold emission of electrons in presence of strong magnetic field, relevant 
for strongly magnetized neutron stars is obtained. It is found that in this scenario, a scalar type 
potential barrier does not allow quantum tunneling through the surface. Whereas, in presence of 
a vector type surface barrier, the probability of electron emission is much larger compared to the 
original Fowler-Nordheim cold emission of electrons. It is found that the relativistic version in 
presence of strong magnetic field does not follow exponential decay. 
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INTRODUCTION 



There are mainly three kinds of electron emission processes from metal surface, they are: (i) thermionic emission, 
(ii) photoelectric emission and (iii) cold emission or field emission. 

Thermionic emission is the external thermal energy induced flow of electrons through a metallic surface barrier 
into vacuum. This occurs because the thermal energy given to the electrons overcome the confining barrier potential. 
Therefore, to have a thermo-electric current, the thermal energy or equivalently the temperature should have a lower 
critical value. This lower limit strongly depends on the work function of the material, which is the minimum energy 
I ■ needed to remove an electron with zero kinetic energy from the material surface. 

The photoelectric effect is a phenomenon in which electrons are emitted from matter as a result of absorption of 
electromagnetic energy in the form of quanta of the radiation incident on their surfaces. Analogous to the thermionic 
emission, in the case of photo-emission, the incident electromagnetic waves or photons must have a lower critical 
'— 1 . frequency, below which there will be no photo-emission. This lower limit also strongly depends on the work function 
1 ' of the material. 

Finally the field emission (also known as cold emission) is an electron emission process induced by strong external 
' electric fields at zero or extremely low temperature. Field emission can happen from solid and liquid surfaces, or from 
I individual atoms. It has been noticed that the field emission from metals occurs in presence of high electric field: the 
t-H ■ gradients are typically higher than 1000 volts per micron and strongly dependent upon the work function. In modern 
science and technology, there are a lot of applications of field emissions, which includes construction of high-resolution 
. and high luminosity electron microscopes and also used in modern nano science and technology 

Further, the thermal emission from metallic surface can be explained by classical Boltzmann statistics, with only 
quantum correction needed is the inclusion of Fermi energy. On the other hand, the photo-emission of electrons 
can not be explained unless the quantum nature of electromagnetic wave is considered, i.e., old quantum mechanics 
, or quantum theory is needed where photon concept of light is introduced. Of course, the concept of work function 
associated with any kind of electron emission process is a quantum mechanical phenomenon. The field emission 
process on the other hand can only be explained by quantum tunneling of electrons. Which has no counter classical 
explanation. However, for general type of surface barrier, this purely quantum mechanical problem can not be solved 
exactly, a semi-classical approximation, known as WKB is needed to get tunneling coefficients. Now because of 
quantum fluctuation, electrons from the sea of conduction electrons (electron gas) always try to tunnel out through 
the metallic surface (surface barrier). However, as soon as an electron comes out, it induces an image charge on the 
metal surface, which pulls it back and does not allow it to move far away in the atomic scale. But if a strong attractive 
electrostatic field is applied near the metallic surface, then depending on the Fermi energy of electrons and the height 
of the surface barrier, it may overcome the effect of image charge and the electrons get liberated. Since the external 
strong electric field is causing such emission and does not depend on the thermal properties of the metal, even the 
metal can be at zero temperature, it is called field emission or cold emission. 

The theory of field emission from bulk metals was first proposed by Fowler and Nordheim in an epoch making 
paper in the proceedings of Royal Society of London in the year 1928 0t3- Fowler-Nordheim tunneling is the wave- 
mechanical tunneling of electrons through a triangular type barrier produced at the surface of an electron conductor 
by applying a very high electric field. Fowler Nordheim tunneling current is found to have important implication in 
various scientific and technological applications, starting from cold cathode emission to VLSI. 

Cold emission or field emission processes not only have significance in the terrestrial laboratories, but is found to be 
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equally important in the electron emission processes from cold and compact stellar objects. For an example, in the case 
of a rotating neutron star, the existing large magnetic field, > 10 12 G for the conventional radio pulsars and > 10 15 G 
for the magnetars, produces a strong electric field at the poles, approximately given by F ~ 2 x 10 8 P _1 Bi2Volt cm -1 , 
parallel to B at the poles 0, here P is the time period of the neutron star in second. At the proximity of polar 
region the potential difference changes almost linearly with distance from the polar cap, which is a region very close 
to the magnetic pole. The repulsive surface barrier in combination with this attractive potential, forms a triangular 
type barrier at the poles. Now, the study of plasma formation in a pulsar magnetosphere is a quite old but still an 
unresolved astrophysical issue. In the formation of magneto-spheric plasma, it is generally assumed that there must 
be an initial high energy electron flux from the magnetized neutron stars. At the poles of a neutron star the emitted 
charged particles flow only along the magnetic field lines. The flow of high energy electrons along the direction of 
magnetic field lines and their penetration through the light cylinder is pictured with the current carrying conductors. 
Naturally, if the conductor is broken near the pulsar surface the entire potential difference will be developed across a 
thin gap, called polar gap. This is of course based on the assumption that above a critical height from the polar gap, 
because of high electrical conductivity of the plasma, the electric field F, parallel with the magnetic field near the 
poles is quenched. Further, a steady acceleration of electrons originating at the polar region of neutron stars, traveling 
along the field lines, will produce magnetically convertible curvature 7-rays. If these curvature 7-ray photons have 
energies > 2mc 2 (with m the electron rest mass and c the velocity of light), then pairs of — e + will be produced 
in enormous amount with very high efficiency near the polar gap. These produced e~ — e + pairs form what is 
known as the magneto-spheric plasma |6l-ll2l| . The cold emission, therefore plays a significant role in magneto-spheric 
plasma formation. In turn, the motion of charged particles in the magnetosphere in presence of strong magnetic field 
causes pulsar emission in the form of synchrotron radiation. Therefore the cold emission process indirectly affects the 
intensity of synchrotron radiation. 

Although the field electron emission is a paradigm example of what we call the quantum mechanical tunneling, 
unfortunately, it is also a paradigm example of the intense mathematical difficulties associated with this tunneling 
process with general type barrier structure. Exactly solvable models with simple type tunneling barrier lead to 
equations 0, y] that underestimates the emission current density by a factor of 1000 or more. If a more realistic type 
barrier model is used by inserting an exact potential barrier in the simplest form of the Schrodinger equation, then a 
complicated mathematical problem with very little physical important arises over the resulting differential equation: 
it is known to be mathematically impossible in principle to solve this equation exactly in terms of the usual functions 
of mathematical physics, or in any simple way. To get even an approximate solution, it is necessary to use special 
approximate methods known in physics as semi-classical or WKB approximation. 

Moreover, to the best of our knowledge, a relativistic version of cold emission model in presence of strong magnetic 
field, relevant for electron emission from the polar region of strongly magnetized neutron stars, even with a simple 
type potential barrier has not been reported in the literature. It has already been mentioned that the emission of 
high energy electrons and thereby formation of magnetosphere is caused by strong electric field at the polar region 
and this strong electric field at the poles is produced by the magnetic field of rotating neutron stars. 

We have organized the article in the following manner: in section-2 we have studied electron emission with a scalar 
type potential barrier. In section-3 we have studied the same thing replacing the scalar potential by a vector type 
potential barrier. In the last section we have given the conclusion and future prospect of this work. 

2. RELATIVISTIC EMISSION MODEL IN PRESENCE OF MAGNETIC FIELD: SCALAR POTENTIAL 

In the original paper by Fowler and Nordheim on cold emission, the conventional non-relativistic quantum mechan- 
ical equations, i.e., the Schrodinger equation was used @. In the relativistic version we shall use the Dirac equation 
and introduce the triangular form of barrier potential at the surface of the system using the conventional form of 
relativistic quantum hadro-dynamics. In the simplified picture of quantum hadro-dynamic model the interaction can 
either be of scalar type or vector type or both. In this model the conventional form of Dirac equation with scalar and 
vector type interactions is given by 13] 

W (idp - gvVJ - (m - g s <f>)} j> = (1) 

where g v and g s are the vector and scalar coupling constants respectively, and <fi are respectively the vector 
and the scalar fields and j^s are the Dirac 7-matrices. In presence of strong neutron star crustal magnetic field, 

electrons can travel almost along the field direction. If the magnetic field strength B > B^ , where B^ is the typical 
strength of magnetic field at which the electron cyclotron quantum exceeds the corresponding rest mass energy, 
then in the relativistic scenario the Landau levels of the electrons will be populated. In this situation, the electron 
momentum in the plane orthogonal to the direction of magnetic field, gets quantized and is given by pj_ = (2veB) 1 / 2 , 
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where v = 0,1,2.... are the Landau quantum numbers. Assuming that the magnetic field is along z-direction the 
corresponding electron momentum component changes continuously from — oo to +00. It can be shown that the 

critical strength Be w 4.43 x 10 13 G. In this situation the Dirac equation reduces to effectively one-dimensional in 
nature. Further, we replace the scalar and vector part of the fields by the non-relativistic potentials at the interface 
between matter and vacuum with triangular in shapes. In the case of vector potential, only the zeroth part is retained 
whereas its spacial part vanishes because of rotational symmetry within the matter. In one-dimension, for the sake 
of convenience we use the form of Dirac equation originally proposed by Dirac, where 7-matrices are replaced by a 
and (5 matrices. The triangular shape potential barrier is represented by the form V(x) = C — eFx, where following 
the notation of Fowler and Nordheim, C is the surface barrier and F (absorbing the electron charge, we express eF 
by F) is the electric field at the polar region. The quantity F acts as the driving force for electron emission. Let us 
first assume that the potential barrier is of scalar type. Then the Dirac equation in one-dimension is given by 

[ap + f3(m + C- Fx)} ijj = Eip (2) 

where E is the energy eigen value. Following the representation [T3] 

i ~o ) and 13 = = ( 1 ) (3) 

and defining m* = m + C as the effective electron mass, then with the spinor notation 



we have 



dx 

After eliminating v from the eqns.(5)-(6), we have 



dv 

— + (to* - Fx)v = Eu (5) 
dx 

'to* — Fx)u = Ev (6) 



^ + Fu + (to* - Fxfu = E 2 u (7) 



dx 

To solve this equation analytically, we make the following transformation of coordinate 



i = F^[^r-x) (8) 



Then we have 



d 2 u 2 / E 2 



This equation may be re-arranged in the following form 

^ + (A~e> = (10) 

which is the well known Hermite differential equation, with A = 1 + 2^, which gives E 2 — 2{v + 1)F and the solution 
is given by 

u = N-^-exp(- l plH v (0 = u I (say) (11) 



F i/2 — * y 2 

where N is the normalization constant and H u {^) is the well known Hermite polynomial of order v (which is not 
necessarily integer) . It is trivial to show that v will also have the same kind of solution. This spinor solution (u or v) 
is for the electron which has already been liberated through the surface into vacuum under the influence of electric 
field F. 
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Now the one-dimensional form of free Dirac equation for the electrons bound within the metal body (in this case 
within the metallic crustal matter), i.e., before emission from the system is given by 

(ap + 13m) ip = Eip (12) 

Then following the same technique as above, we have 

— + w lu = (13) 

where Wk = (E 2 — m 2 ) 1 / 2 (exactly same type of equation can be obtained for v). Now following the notation of Fowler 
and Nordheim, we express the solution for free electrons within the system, confined by the surface barrier V(x), in 
the form 

[aexp(iwkx) + a' cxp(— iwkx)} = ujj (say) (14) 
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where a is the probability amplitude for electrons moving along the positive direction of x, whereas a' is the corre- 
sponding quantity for left moving waves. Assuming the interface between the neutron star matter (which in this case 
is the dense metallic iron crystal of neutron star crustal matter) and vacuum is at x = 0, the wave function and their 
derivatives must be continuous at x = 0, i.e., 

Ujj(0) = uj(0) and m^-(O) = u'j(0) (15) 

Using the relation H' v (£) = 2vH v -\{£), we have 

a + a' = Nwl /2 j^ cxp (-|) H^ ) (16) 

and 

iwl /2 (a - a') = NEexp (-^j [S H V (£ ) - 2uH v - 1 (£ )] (17) 

where = £,( x = 0) = m*/F 1 ^ 2 . These two conditions (eqns.(16)-(17)) may be rearranged in the following form 

a + a' = X and a - a' = iY (18) 
where X and Y are two real quantities. Hence it is straight forward to verify that the transmission coefficient, given 

by 

Tr=1 -w> (19) 

vanishes exactly. In this section we therefore conclude that if the barrier in combination with the external electro- 
static driving force behaves like a scalar potential, electrons can not tunnel through the surface barrier whatever be 
their kinetic energies and the strength of external electric field. Identical result will be obtained if v is used instead of 
u. This result has some analogy with the confinement of quarks inside hadronic bag with scalar type linear potential. 



3. RELATIVISTIC EMISSION MODEL IN PRESENCE OF MAGNETIC FIELD: VECTOR POTENTIAL 

In presence of a vector type potential barrier at the interface with the same kind of triangular shape as mentioned 
in the previous section, we have the conventional form of Dirac equation 

(ap + Pm)tp = {E-C + Fx)ip (20) 

Decomposing the spinor ip into u and v as before and using the same representation for a and /3, we have 

— +mu = {E-C + Fx)v (21) 

fill 

+ mv = (E-C + Fx)u (22) 
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Eliminating v from these two equations (eqns.(21)-(22)) we have 

d 2 u [du/dx + mu 



m 2 u ={E-C + Fx) 2 u (23) 



dx 2 ' E-C + Fx 
Substituting £ = — (C — E)/F + x, we have 

^ (mC + ~ F ^ 4) " = ° (24) 

To obtain an analytical solution for u, let us consider the non-linear transformation 77 = £ 2 . Then the above equation 
becomes 



d 2 u 1 / /2\ 7 1/9 2r7m 2 , 

»? 7 + -7T- U = (25) 



1 / /ox 1/2 

(. -K 

drj 2 Ari 2 \ "~ \F J ' ' F 

To obtain a closed analytical solution, we neglect for the sake of simplicity the rj 1 ^ 2 term. Then we have 

d 2 u / 2x 



, 1 - « = (26) 
dr\ \ T] ) 

where x = rn 2 /AF and 77 > 0. This is the standard differential equation for the well-known Coulomb function with 
zero angular momentum (l5j . The solution is given by 

u = —5- exp(2z<5 ) exp[i(r7 - xln 2tj)} = m (say) (27) 

2rj 

where <5q = argr(l + i\). 



The exact solution is quite complicated. A semi-compact solution may be obtained using Mapel package in terms 
of HeunB function- the Heun Biconfluent function. This function is related to Kummer function or Hypergeometric 
function. The solution is given by: 

u(rf) = C\ exp( i?7)HeunB(2, 0, -2i X , (4 + Ai)p, (1 - i)r] 1/2 )r) + C 2 exp^HeunB^, 0, -2i X , (4 + U)p, (1 - iW /2 )v 

exp(— 2irj) 



(HcunB(2, 0, -2i X , (4 + U)p, (1 - i)r/ 1 /2))2 T? : 



where 



/ 1 \ V2 
m / 2 x 



9 4 \F 



To obtain tunneling co-efficient , we again consider the continuity of u and its derivative at the interface (i.e., at 
x = 0). Which means 

ujj(0)=u/(0) and u' n (0) = u'j(0) (28) 

Hence we get 



(a + a ) = w k < — (29) 



where 

X = Co exp(2i(5 ) exp[i(?7o - x m 2 ?7o)] (30) 
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is the phase factor, and 



/a X . 
(°-°) = 1/2 i 
V w k 



d-^) + i- 



Where, as before 770 = T)(x — 0). Hence, from the definition of tunneling coefficient, we have 

\a'\ 2 



l + Vo 



f +% /2 (l- 



(31) 

(32) 
(33) 



Substituting 770 = (C — E) 2 /2F, we finally get 



2 5/2_pl/2 



T B = 



(C-£) 3 (l 



8F 2 + (C - E) 2 uokF 1 / 2 + 2V2(C - E) (1 - 



2(C-B) i! 



(34) 



Obviously, in this case, unlike the scalar type surface barrier, the transmission co-efficient is non-zero. Further, 
unlike Fowler Nordheim non-relativistic model, here the tunneling coefficient dose not fall exponentially. In the above 
expression for tunneling coefficient, as before C = (i + Wf and u>k — (E 2 — to 2 ) 1 / 2 , with m„ = (m 2 + ZveB) 1 / 2 

and v — 0, 1,2 are the Landau quantum numbers for the electrons. In presence of magnetic field, it has been 

observed that the work function increases with the increasing field strength. A numerically fitted form for magnetic 
field dependence is given by W f = W C (B/Bi e) ) - 5 in eV, where W c is a real constant parameter rs 82.93 [U EH. 
One can easily verify that the tunneling coefficient will vanish for F = and extremely large values for electron work 
function. 



4. RELATIVISTIC FORM OF COLD EMISSION CURRENT FOR B ^ 

To obtain electron cold emission current in the relativistic scenario in presence of strong magnetic field, we consider 
the motion of relativistic electrons along the direction of magnetic field and the system is assumed to be at zero 
temperature. It is well known that in presence of strong quantizing magnetic field > 10 15 G, which is the recently 
observed surface values of some strange type neutron stars, known as magnetars 16], the motion of the electrons at 
the crustal region of a neutron star is effectively one dimensional in nature; motion is restricted along the direction 
of flux lines. In this case the density of states is given by [l9| 



d 3 P _ eB ^ 
(2tt) 3 4tt 2 ^ [ 



Svo)dp z 



(35) 



In this expression we have taken Planck constant h = 1, velocity of light c = 1 and Boltzmann constant fc# = 1; 
known as natural unit. Here v is the Landau quantum number and the factor 2 — S u q takes care of singly degenerate 
zeroth Landau level and doubly degenerate other levels with v ^ 0. The upper limit of Landau quantum number is 
infinity at finite temperature, whereas at zero temperature it is given by 



2eB 



(36) 



here /i is the electron chemical potential and [ ] indicates that v max is an integer but less than the actual number. 
Now the number of electrons incident normally on unit area at the interface between matter and vacuum is given by 



*=SB 2 -M/~f— 

^ V=Q Jpo exp 



dp z 



(37) 



where E v — (p 2 + to 2 + 2^eB) 1/ ' 2 is the electron energy at ^th Landau level and po is the lower limit of z component 
of electron momentum. Unlike thermionic emission and photo-emission, the cold emission is a purely quantum 
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mechanical phenomenon; electrons tunnel through the surface barrier and its emission in free space is driven by the 
external electric field. Since the emission process is controlled by the tunneling probability of electrons through surface 
barrier, the actual expression for cold emission current will therefore be 



2f> 00 rc 



Pz 



dp z 



E, 



v exp 



T 



(38) 



In the case of neutron star, since the temperature is low enough compared to the chemical potential of electrons, we 
put T = in the Fermi distribution in eqn.(38), then we have after replacing the integration variables from p z to E v 
and omitting the index v (keeping in mind that E = E{p z ,i>, B)) in electron energy 



R 



e 2 B 
4^ 



I' max p/j, 

^ (2 - M / T R dE 



(39) 



We have obtained cold emission current by evaluating numerically energy integral and the sum over Landau quantum 
number. However, we have noticed that at the crust of magnetars the Landau quantum number becomes exactly 
zero, i.e., the electrons occupy only their zeroth Landau level. In fig.(l) we have shown the variation of cold current 
with electric field and magnetic field by the solid curve and the dashed curve respectively. Both the curves indicate 
the initial increase of cold current. Since the work function increases with increasing magnetic field and electric field 
is generated by the rotating magnetic field of neutron stars, the cold current decreases abruptly after some magnetic 
field / electric field strength. 



5. CONCLUSIONS 



In this article we have studied cold emission of electrons from strongly magnetized neutron star polar region using 
rclativistic quantum mechanics. To get an analytical solution, we have considered triangular type potential barrier as 
was used by Fowler and Nordheim. We have noticed that a scalar type potential barrier does not allow electrons to 
tunnel out the interface, whatever be the strength of external electric field. With the vector type potential barrier, we 
have obtained an analytical solution for the tunneling coefficient. Unlike the original work of Fowler and Nordheim 
the tunneling coefficient does not follow any exponential law. However, the variation of cold current with the electric 
field strength can be fitted numerically using a x 2 -minimization program. The numerically fitted functional form is 
given by 

R = 0.26F 2 4 /2 exp(-9.8 x F14) (40) 

where F u = 10~ 14 F and F 2 a = 10~ 24 F. The variation shows an exponential decay form is because of approximate 
linear dependence of electric field with the magnetic field strength of the rotatiting neutron star. We have obtained 
an analytical expression for the cold emission current which changes with the strength of magnetic field also. We 
found that the cold current initially increases almost linearly with the external electric field and finally when the work 
function becomes large enough because of high magnetic field strength, it drops abruptly and finally goes to zero. 
Since the emission current becomes extremely small in presence of ultra-strong magnetic field, the magnetosphere of 
the magnetars will be of very low in charge density. As a consequence there will be very weak synchrotron emission 
in the radio wave band. 
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FIG. 1: Variation of cold emission current with neutron star magnetic field (dashed curve) expressed in terms of critical field 
for electron and electric field in Volt/cm (solid curve) at the neutron star pole. 
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